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The measuring process is studied, where a macroscopic number N of particles in the detector 
interact with the object. When N — » oo, the fluctuation of the object between different eigen-states 
of the operator O to be measured is suppressed, frozen to one and the same state while the detector 
is on. During this period, the stationary phase accompanying the macrovariable is established to 
have a one to one correspondence with the eigen- value of O. A model is studied which produces the 
ideal result when N — > oo and the correction terms are calculated in powers of 1/JV. It is identical 
to the expansion including the fluctuation of the object successively. 
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I. INTRODUCTION 



The observational problem in quantum mechanics has 
a long history of debates [3, S, H 0, B IE 0, H ■ In par- 
ticular, the role of the docoherence due to environmental 
influence has been widely discussed [E EE EE EE EH- 
The dynamical reduction model has actually been con- 
structed and many recent researches are focused on this 
subject [3, EE [3 EE 13 • Irrespectively of the mecha- 
nism that leads to the reduction, we need the statistical 
treatment of the data, which is done by using the wave 
function following the rule of the ordinary quantum me- 
chanics. When we apply the quantum mechanics to the 
detector system, composed of a number of microscopic 
particles, two requirements have to be fulfilled; (1) any 
detector variable X(t), the pointer position for example, 
should show a non-fluctuating behavior of classical type 
as a function of the time. (2) different eigen-states of 
object operator O to be measured have to be mapped 
onto different values of X(t). 

We stress here that both (1) and (2) are realized by 
the stationary phase accompanying any macroscopic vari- 
ables [III] (The precise definition of the macrovariable is 
given in Sec lII A| ). Consider (1) in the path-integral form 
[20l | . Out of many fluctuating paths, the stationary phase 
selects one particular smooth path denoted as X st (t). 
Then the absolute square of the wave function of X inte- 
grated by other degrees has a non-diffusive peak, equal to 
the density of a classical point-like particle 8(X — X st (t)). 
This is what we observe as a macroscopic object in the 
ordinary life. As for (2), we recall that the measuring de- 
vice is prepared in such a way that O interacts, directly 
or indirectly, with a large number N of particles in the 
detector. So the Hamiltonian Hi describing such an in- 
teraction may be a function of O, coordinates Xj's and 
the momenta ^'s of all these particles and is O(N). Since 
the Hamiltonian of the object Hq itself is 0(1), it can be 
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neglected compared with Hi. Then, for N oo, each 
eigen-state of O is frozen in the same state as long as the 
detector is on. Moreover, precisely during this period, the 
detector variable X(t) changes its value depending on the 
eigen-values of O. 



The desired mapping is realized in this way. Note that 
if the object interacts with a finite number of particles, 
the process is simply a quantum mechanical one, nothing 
to do with the measurement. 



The above observation suggests the 1/JV expansion 
scheme for large N, which incoorporates successively the 
fluctuation of the object connecting different eigen-states 
of O. It is the purpose of this paper to show that this is 
indeed the case by adopting a simple model of separable 
X. It is solvable by 1/JV expansion and we calculate sev- 
eral lower order terms explicitly. The stationary phase 
was applied to the macrovariable in [l9j in the lowest 
order but above observations were lacking. 



In Sec|lU crucial points of the stationary phase are pre- 
sented on which subsequent discussions are based. These 
are not stated in [l9j]. Following Sec (Til a model is de- 
fined in Sec lIIIi and the signal function is introduced, 
which agrees with the density of the classical point-like 
particle for TV — * oo. Higher order corrections are stud- 
ied in Sec lIVI They are given by the expansion in (frac- 
tional) powers of 1/JV and the result is summarized in 
(|85p below. This type of calculation is surely required 
since in the actual detector N is finite, although very 
large, and the correction terms might be tested experi- 
mentally. An attempt at the numerical estimate of the 
correction terms is presented. After the discussions in 
Sec[V] general non-separable case is discussed in Ap- 
pendix [SJ In Appendixes IBICID1 some of calculational 
details are shown and in [El the order estimation of ne- 
glected terms is given. 
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II. STATIONARY PHASE AND 
MACRO VARIABLE 

A. Macroscopic system and macovariable 

Before constructing the model of measurement, the 
properties of macro and microvariables have to be elu- 
cidated, since the quantum mechanical detection pro- 
cess consists of an interplay between the two. Con- 
sider a macroscopic system, which contains a large num- 
ber of microscopic particles with coordinates Xi (i — 
1, 2, • • • , N), all of which are assumed to have the same 
mass m. To make formulas simple, we work below in 
one dimension, extension to three dimensional case be- 
ing straightforward. Now an extensive quantity grows 
up with the size of the system. As has been introduced 
in [l9j in the case of field theory, a macrovarible of N 
particle system is an intensive quantity defined by di- 
viding an extensive quantity by TV. The center of mass, 
X = X^ili X i/N is the simplest example. 

To make above definition more precise in quantum 
case, let us recall that the quantum mechanical process 
is described by the path-integral form given by 




where S = j dtL is the action functional of the system. 
The Lagrangian L = L([xi], [xi]) is an extensive quantity, 
so it is O(N). (([xi] represents Xi's collectively.) For this 
statement to hold, the system has to be thermodydnam- 
ically normal, which holds when the interation among 
microdegrees is of short range and the particle density is 
finite over all space. (The statistical factor 1/Nl has to 
be inserted properly for the system of identical particles.) 

The macroscopic system behaves classically as a whole, 
while the system at the same time contains atoms and 
molecules, which is described quantum mechanically. 
The stationary phase accompanying the macrovariable 
just realizes such a picture. The reason is simple; if we 
change one X4 as Xi — > x% + a by a finite amount a, then 
the change of the phase S[xi] is of the order a. On the 
other hand, suppose the macrovariable X changes by a 
finite amount a. It means that macroscopic number N of 
Xj's are displaced by the order of a. Then S[xi] in general 
shifts by O(Na). Now we integrate over all values x^s in 
([]]), so when N — > 00 only that point of S[xi] contributes 
which is stationary when X varies. Since this holds for 
any time slice of the path- integral ((T|) , one gets a smooth 
stationary trajectory of X. It describes the motion of a 
macroacopic body as a whole. On the other hand, no 
stationary point exists for each Xi, so every Xi remains 
as a fluctuating quantum variable. Consider instead the 
limit H — ► 0. The change in S[xi]/H under Xi — > Xi + a 
is O(afh) for each i, so every Xi is determined by the 
stationary equation, i.e. the Newtonian equation. It is 
not the "classical limit" as seen in the ordinary life. 



B. Particle picture by the stationary phase 

The above statement is seen clearly if we take the sep- 
arable case of the center of mass for X. (See Appendix lAl 
for non-separable case.) We write Xi = X + x\, where x' t 
is the coordinate measured from the center of mass and 
satisfies J2i x i = 0- Now |T]) becomes 

J[dX}e iS ^/ h J {[[dx' l }8(Y, x 'iy S2lX ' i]/h - ( 2 ) 

i— 1 i 

The second factor describes the microscopic quantum 
phenomena and can be neglected, or integrated out, in 
the measurement thory, since we are interested only in 
X. The first factor accompanies the stationay phase 
since, as discussed above, Si [X] is proportional to N; 
Si[X] = Ns±[X]. For large N, the X integration is dom- 
inated by the solution of functional stationary equation 
5s\[X]/8X{t) = 0. The fluctuation of X is suppressed 
by the phase cancellation and a single smooth trajectory 
X st (t) is selected by the constructive phase coherence 
among paths near X st (t). Although the fluctuation of 
X defines the wave function *&(X, t) for finite N, once 
N becomes infinity, X reduces to the variable of a point- 
like particle. These statements are based on the following 
well-known formula. Let f(X) be a function of X having 
the stationary point at X$, then we have 

^J^=^^Z 5( X-X„). (3) 

Consider here the Feynman kernel K(X, T; Y, 0), which 
connects the wave functions at different times; 

V(X,T) = J K{X,T;Y,0)^{Y,0)dY. (4) 

Applying ([3]) at every time slice from t = to t = T, 
the kernel is seen to contain a factor 5(X — X st (Y,T)), 
where X st (Y, T) passes Y at T — 0. (Initial velocity 
depends on the form of ^(X, 0).) Thus, each point X 
on the wave function ^>(X, 0) is just transported along 
X st (Y, T), as opposed to the Huygens picture of wave 
mechanics. Our ty(X, T) here represents the mixed state. 
If we choose |*(F,0)| 2 = 5(Y - X x ), then \<$(X,T)\ 2 = 
5 (X — A" st (A"i, T)); the wave function has a non-diffusive 
sharp peak, representing the density of a point-like par- 
ticle. Explicit examples appear later. 

To discuss other local densities, let us discretize the 
time with the interval At. Then the fluctuating momen- 
tum operator [H P = (h/i)d/dX{t) = Nm(X(t + At) - 
X{t))/ At becomes the classical expression NmX(t) eval- 
uated along X 8t (t) as N — > 00. The quantum mechanical 
expression of the momentum or the energy e.t.c. reduces 
to the corresponding classical density; 

V*(X,t){ 1, P, P 2 /2Nm )&(X,t) — ► 
( 1, NmX(t), NmX 2 /2 )8(X - X st (t)). 
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C. Measurement by macrovariable 

Suppose we measure the object operator by the in- 
teraction Hamiltonian Hj(0, [xi], \pi])- Then the total 
Hamiltonian is the sum of three terms; 

H = H + HnQxt], fa]) + H T (0, [x t ], fa]), (5) 

where Hq is the object Hamiltonian, Hp that of the 
detector. As stated in the Introduction, O interacts with 
many particles in the detector and X is so chosen that it 
includes (almost) all of them; 

JV JV 

H! = J2HO,Xi,Pi), X = Y,9{xuPi)/N ( 6 ) 

8=1 1=1 

with some functions h and g. Here and hereafter, N is 
the number of particles involved in ©. Hn in (JSJ) is 
the Hamiltonian of these particles and is O(N) together 
with Hj. Arguments of Sees III Al III Bl can be applied to 
X thus defined. 

One comment here. If one can find any parameter a 
which produces the stationary phase for some variable rj, 
then a and rj can be used in place of N and X. Then, 
the macrocity will not be required for the measurement. 



D. Freezing the object state 

Now in general, O does not commute with Ho, 

[H o ,O]?0. (7) 

So, the eigenstates of O always fluctuates among differ- 
ent states. To measure such a fluctuating object, we have 
to suppress the size of the fluctuation somehow. This is 
done by taking Hj = O(N); since Ho = 0(1), we can 
neglect Ho for N — > oo and adopt H = Hq + Hj as 
the Hamiltonian. Thus, taking the representation which 
diagonalizes O, 0\a >= X a \a >, each eigen-state \a > 
develops by H D + Hi(X a , [x^, fa]) = H D + H 1>a and re- 
mains in the same state \a > as long as Hj ^ 0, i.e. the 
detector is on. Precisely during this period of Hi ^ 0, 
the macrovariable can change its value since it evolves by 
Hd + Hi, a , which produces different stationary phases 
for differnt A a 's since Hi is O(N). Thus we get different 
values of X depending on the microscopic state of the 
object. It is realized independently of the detailed forms 
of the functions h or g above, as long as both include 
large number of particles which interact with O. This is 
the amplification mechanism of the detection process in 
terms of the stationary phase. Higher order correction 
terms for large N are given systematically through the 
expantion in powers of the off-diagonal elements of the 
object Hamiltonian Ho- Below, this is done explicitly 
based on the model Hamiltonian. 



III. THE MODEL 
A. The Hamiltonian 

When we construct the model, above functions 
h(0,Xi,pi) and g(xi,pi) has to be fixed. The simplest 
case is h = —fOxi and g = Xi (/ is the coupling 
strength.); then 

N 

H z = -J2 fOxi = -fNOX. (8) 

i=i 

Here X is the center of mass of particles that interact 
with O. In the realistic detector, the photomultiplier 
for example, complicated processes may happen. An ob- 
ject interacts with an atom in the detector (via exchang- 
ing a photon), ionizing an electron. It is accelerated by 
the electric field applied in certain direction, which inter- 
acts with another atom, ionizing a second electron and 
so force, until we have a macroscopic number of elec- 
torns, giving a signal as the current. Or a high energy 
object interacts with many atoms in the direction on its 
momentum, along which the track of ionized electrons is 
detected. Above Hi simulates these processes by a di- 
rect interaction of O and N particles, which are in the 
direction of applied electric field or in the direction of the 
object momentum. So the problem can be simulated by 
one dimension, with X taken to be the center of mass of 
Xi in that direction. Independently of the detailed form of 
the interaction, the essential point is that, although each 
electron receives a microscopic amount of energy, the sum 
of them is O(N), which affects the stationary phase of X. 
Since H D is P 2 /2Nm = NmX 2 /2 (P = NmX is the to- 
tal momentum) plus terms independent of X, the total 
Hamiltonian, the object plus detector, of our model is 
defined by 

H = H + NmX 2 /2- NfOX = H + H N , (9) 
H N = NmX 2 /2 - NfOX. (10) 

We do not write x\ part hereafter, since the dependence 
on X and x\ is factorized as ©. 



B. The initial wave function 

Below the eigen- values A Q of O are assumed to be dis- 
crete and non-degenerate; A a ^ Xb if a ^ b. Writing the 
eigenstate of X as \X > with continuous eigen- value X, 
the complete set of states of our model Hamiltonian is 
given by \a > \X >= \a,X >. We also use the com- 
plete set spanned by \x,X >. The eigen-function is then 
4> a {x) =< x\a >. 

Let the detector be switched on at t — 0, and the state 
vector at t = is written as \^ >= \<f> > |<& >, which is 
the product of the object \<j> > and of the detector |$ >. 
Expanding as \(f> >= C a \a > into complete sets, the 
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initial wave function is given by 

V(x,X,t = 0) =< x,X\V >=< x\<f> >< X\<f> >, 

< X\(j) >= J2 Ca < X \ a >= ^2 C a4>a{x), (11) 

a a 

The initial wave function of the detector *&(X) is assumed 
to have a peak at some value of X , with the precision A. 
To be explicit, we adopt a Gaussian type; 

< X\& >= (ttA 2 )" 1 / 4 exp(-X 2 /2A 2 ). (12) 

^(X)] 2 becomes S(X) when A — > 0. The center position 
is at X — and the initial velocity X — (h/Nm\)d/dX is 
also zero when N — ► oo. (Non-zero velocity v is obtained 
by multiplying exp(iNmvX/h) to ([T2"]) .) 

When N — ► oo, the center of the peak traces a clas- 
sical trajectory determined by Hn keeping the width A 
constant. Here <X\^f> | 2 represents the density matrix 
of a mixed state, with A repesenting the classical un- 
certainty. If N is large but finite, X fluctuates and the 
diffusion process comes in. The numerical consideration 
is given at the end of Sec lHI El where we will see that, 
up to the order we are considering, the influence of the 
diffusion is negligible in the detection process and the 
dominant effect comes from the fluctuation of the object 
while the detector is on. 

Now the macroscopic limit is N — > oo, with other quan- 
tities kept fixed. But in order to avoid the classical un- 
certainty and keep various formulas simple, we take the 
limit A — > of the coefficients of limitting expression 
obtained by N — > oo. Such a limit, first N — > oo, then 
A — > 0, is denoted as 



C. Time evolution and the expansion by the off 
diagonal elements of Ho 

The total wave function develops in time as 

X, T) =< x, X\ exp(-LHT/7i)|* > 
=< x,X\exp(-iHT/K)^C a \a > |$ > . (13) 

a 

When we expand in powers of Ho, we first sum up the 
diagonal term (Ho)aa exactly in every order of expan- 
sion. Then the expansion becomes the one in terms of 
the power of off-diagonal elements (Ho) a b, a ^ b. To 
achieve this, we use the well-known formula; 

U{T) = exp(-iHT/h) = U N (T)U (T), 
Un(T) = exp(—iH N T/h), 

U {T) = Texp(^-i^ dsH {s)/h\, (14) 
H (s) = u]f(a)HoU N (a). (15) 



In (fl4| . T implies the time ordering operation. To get 
the desired expansion, consider 

H (s)\cl >= tf N (s)H U N {s)\a > 
= U^ a (s)H \a>U N<a (s) 

= U N, a ( s )( H o)aa\a> U N: a{s) + off diagonal part 
= (Ho)aa\a > +off diagonal part. (16) 

Here we have introduced 

U N ,a( s ) = exp(-iH N . a s/H), (17) 
H N ^ a - NmX 2 /2 - Nf a X, f a = fX a . (18) 

Thus, when we sum up the diagonal elements of Ho, Ho 
can be treated as a c-number and the diagonal parts are 
summed up into the phase. Thus we can use the formula 

U N (T)U (T)\a >= U N , a (T) exp(i0 a T)|a >, 

0a = -{H )aa/h. (19) 

Now we concentrate on the defining equation of the 
T-product (HU). It is an infinite product of the term 
exp(— ii?o(s) As/h) in the infinitesimal time interval As. 
When it is evaluated by off-diagonal elements [i?o] nd , 
we write exp(—i7?o(s) As/ft ^ 1 + (— iHo{s)As/h. Using 
cqs. p^)) and (fl^)) . the expansion thus obtained becomes 

oo 

< b, X\U N (T)Uo(T)\a, Y>=Y, u£\x, Y), (20) 

fe=0 

where involves A:-th power of [Ho] nd - The results of 
first few terms are shown below; ( [Uo(T)] d contains the 
diagonal parts only, while [i?o] nd off-diagonal parts. ) 

U^(X,Y)=<b,X\U N (T)\a,Y > 

= 5 ab cxp(i6 a T) < X\U N , a (T)\Y >, (21) 

c/ fc ( Q 1) (x,r)=<&,x|c/ Ar (T)(-i/ft) 

x f ds[U (T - S )} d [H (s)} nd [Uo(s)} d \a,Y > 
Jo 

= {-i/K)(H )ba [ dsexp{i(9 b {T-s) + e a s)} 
Jo 

x < X\U N , b (T - s)U N . a (s)\Y > . (22) 

In (|22|) and in what follows, ([Ho] nd )ba is written simply 
as (Ho)ba for the notational simplicity so a ^ b is im- 
plied. The factor S ab in (f2"Tj) implies that the object does 
not fluctuate in the lowest order while the detector is on. 
By a similar manipulation, 

Ug\X,Y) = (-i/h) 2 J2(Ho)bc(Ho)ca 

C 

x [ ds' [ ds exp{i((9 fc (T - s') + 9 c {s' - s)) + 6 a s] 
Jo Jo 

x < X\U N<b (T - s')U N , c {s' - s)U N , a (s)\Y > . (23) 
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The wave function has the corresponding expansion 
*(x, X, T) =< x, X\U N {T)U (T) ^ C a \a, $ > 

a 

/OO 
dY^2< x\b > C a J2u£\x,Y) < Y\* > 

6,0 

oo oo 

= E < ^ > E * (fe) ( 5 ' X ' T ) = E * (fc) (*> x > T )- 

& fc=0 k=0 

The wave function in the \a > representation 

oo 

V(b,X,T) =< b,X\U(T)\V >=^2^ {k) (b,X,T) (24) 

fc=0 

has been introduced and each ty( k > has fc-th power of 
[H ] nd . In case [H ,0] = 0, (# ) 6a (A - A 6 ) = 
follows, so one gets [H ] nd = 0. Thus, in (2© or p3]) 
only the lowest term with fc = is non- vanishing. 

D. The signal function 

In the real experiment, the object is not actually ob- 
served, so let us define the signal function by integrating 
(summipng up) |\1/| 2 by x (b); 

J{X,T)= f dx\*(x,X,T)\ 2 = J2\*(t>,X,T)\ 2 (25) 
J b 

oo 

= £jW(X,T), (26) 
j(°\X,T)=Y,\y {0) (b,X,T)\ 2 , (27) 

b 

jW(X,T) = ^*(°)*(&,X,T)#«(6,X,T) + c.c., (28) 

b 

(X, T) = j} 2) (X, T) + 4 2) (X, T) , 
j[ 2) (X,T) =^#)*(6,X,T)*( 1 )(6,J,T), (29) 

6 

J 2 (2) (X,T) =J]*^*(&,X,T)*( 2 )(6,X,T) + c.c.. (30) 



Here J^(X,T) is of the order of ([iJ ] nd ) fc . Let us cal- 
culate J( '^ 2 )(X,T) and l^ * 1 ' 2 ) successively. 



E. * (0) (b,X,T) 
The lowest term is calculated by ([T9"]>. ([21]) and (jUJ as 

*C°J(6,X,T) = < b,X\U N {T)\a,$ > 

a 

= J dY exp(i0 b T)C b < X\U N>b (T)\Y >< V|$ > . 



Now we insert (fT2l . and use the following result of the 
Feynman kernel [2Q| for the Hamiltonian H^.a of (fl8| : 



< X|C/^ (T)|Y >= 



e Q = 



27rifiT 

m(X-y-e a (T)) 2 
2T 



= /aT 2 /2m, 



(31) 



(32) 



£ a (T) is the classical change of X during T in the presence 
of the constant force f a , with the initial condition X — 
X = 0. 8 a is the classical action for this motion. By ((3]), 

lim < X\U N , a (T)\Y >= 5(X-Y- £ a (T))e i§ °. (33) 

N— >oo 

Apart from the phase, the whole wave function develops 
as a parallel transport; 

$W(x,X,T) = e i§ «^°\X - £ a {T),0). (34) 

For the discussions below, the ^-integration is done for 
general N. After a simple Gaussian integral, we get 



^°\x,X,T) = C b e ie " T 

N 2 m 2 {X-UT)) 2 



1 



1/4 



Nm 



x exp 
D = 



T 2 h 2 2D 
1 iNm 
A 2 ~ ~KT' 



\Nm(X-j b (T)f 
2HT 



Taking the absolute square, 



|^(A,T))| 2 = |C h | 2 



Here we have written 
DD* N2m2 



1 



7rA 2 p 



exp 



A 2 p 



T 2 h 2 



h 2 T 2 



P 



(35) 



(36) 



(37) 



The peak of |\I>( ) (b, X,T)\ 2 traces the classical trajectory 
X = £ b (T) and for large N,p = 1+0(1/N 2 ), so the effect 
of the broadening of the width due to the fluctuation of 
X, is 0(1/ N 2 ) so setting p = 1 is may be allowed. See 
blow for the numerical study. 

Now, in order to map microscopically different chan- 
nels | a | b > into a macroscopically distinguishable 
state, A has to be suffiently small compared with the dis- 
tance of the different peaks; | £ a (T) - £ 6 (T)| » A. This 
is the requirement for the detector, which is assumed to 
be the case. To get the numerical value of T required for 
producing a signal, we estimate 



\UT)-Zb(T)\ 



fgT 2 
2m 



» A. 



Thus we get T^>> ^2raA//„. Let a be the atomic scale 
length, then f a a = /A a a ~ fOa is of the atomic energy 
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size (inonization energy, for instance) due to the interac- 
tion between the object and one particle in the detector. 
If we set rather arbitrarily f a a =1 cV=1.6 x 10 _12 erg, 
A = 10~ 3 cm and take a = 10~ 7 cm, then for the case of 
the electron (m =9.1xl0~ 28 g), 



l2mAc 
f a a 



10 



-12 



sec. 



(38) 



This is quite a small number, which does not change 
much for the proton (m = 1.7xl0 _24 g) and for somewhat 
larger or samller A. 

Next, we try to estimate the magnitude of N, for 
which the diffusion process during the measuement time 
T can be neglected. By (131))) , p — 1 « 1, or equiva- 
lently N » hT/(A 2 m) has to be satisfied. This fol- 
lows also from the uncertainty relation. If we adopt 
A = (1(T 3 , l(r 4 )cm, then we get N » (10 6 , 10 8 )T for 
the electron, JV» (10 3 , 10 5 )T for the proton. (T is mea- 
sured in sec.) Since T given in (|38p is quite samll, we 
conclude that the diffusion effect in the detection pro- 
cess is totally negligible. We set p = 1 hereafter. When 
we include higher orders of [i?o] i non-trivial constraint 
on N will emerge, see Sec II V Hi 



F. The signal function in zeroth order 

j(°)(X,T) is given in ([27]), In the limit =>, we get the 
ideal situation in the measurement. Denoting by — ► the 
time evolution after the detector is switched on, we get 

j(°)(A,0) = ^|C* fc | 2 <5(A) = 5(A) (39) 

b 

— J^(X,T) = J2\Cb\ 2 S(X-UT)) (40) 

b 

Note that the signal function becomes the classical den- 
sity of a point particle moving along £h(T). By ideal, 
we mean that the above result is in conformity with the 
usual the quantum mechanical rule; integrating by X in 
the samll region R b = (&(T) - S, &(T) +5),(5>0), 



jW(X,t)dX 



\C b \ 2 . 



(41) 



Rb 



We say that it is the probability for X to stay in R b , 
which in turn implies that the probability of the object 
to be in the state \b > is | Ct, | 2 , since the mapping \b >^ 
£i>(T) is one to one by the stationary phase mechanism. 

Written by the wave function symbolically, the ideal 
measuring process is expressed as 

< x\a > \/S(X) — ► 

a 

*(a;, X, T) = C a < x\a > e^ T+§ ^y/S{X - &(T)). 

a 

Thus the object stays in the same state. (Suare-root of 
the delta-function is ill-defined so we need some regural- 
ization.) 



IV. HIGHER ORDERS 

A. * (1) (6,X,T) 



As is given in (J22J), < X\U Ntb (T - s)U N , a {s)\Y > has 
to be evaluated. This is the evolution kernel defined by 
the Hamiltonian H(t) = NmX 2 /2 - Nf(t)X, where 



/(*) = 



fa j 

fb ; 



for < t < s, 
for s <t<T. 



(42) 



Now we apply the formula for this process [2C 



<Y\U N , b (T-s)U Nta (s)\X >= 
m(Y -X) 2 



Sba — 



2T 



Y 
1 



dtf(t)t - 



X 



2^T exp [n NSba 



dtf(t)(T-t) 



Using 



T jo 

dt [ dt'f(t)f(t')(T-t)t'. 
we get after a straightforward calculation, 



ba 



2T 



Y-X-Z ba (T,s] 



(43) 



+ XQ ba (T, s) + P ba (T, s), 
Q ba (T,s) = f a s + f b (T-s), (44) 



P ba (T,s) 



f 2 

J a 

6m 



s 2 (3T- 2s) 



f 2 
bm 



fafb 

2m 



s(T-s) 2 , (45) 



6a(T, s) = A (2 T S - ,s 2 ) + A(( T - s) 2 . (46) 
2m 2m 

Here, t; ba (T, s) is the classica change of X during T un- 
der the force f(t). Note that & a (T,0) = £ b (T) and 
£ ba (T,T) — £ a (T), in comformity with the fact that 
s is the time Ho acted, making the transition from 
| a > to \b >. One can confirm again that P ba (T,s) + 
XQba(T,s) coincides with the classical action integral 
along £ba(T,s). Using (fT2^) . (f2"4"]l and applying the oper- 
ation J dYJ2 a Ca < Y\<& >, we finally obtain 



* (1) (fo,A,T) 



(47) 



1 



ttA 2 

T 



1/4 



Nm 



^2(-i/h)(H )baC a 



2mhT\l D 

a 

[ ds x expi{0 a s + 9 b (T - s)} 
Jo 

N 2 m 2 R ba (X,T,s) 2 NmiR ba {X,T,s) 2 



x exp 



2h 2 T 2 D 



2KT 



+ -N{XQ ba (T,s) + P ba (T,s)} 
n 

R ba (X,T, S )=X-^ ba (T,s), 
Note that R ba {X, T, 0) = X - ft(T). 



(48) 
(49) 
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B. Signal function J (1) (X,T) 



By ([28]). ([35| and gH]), one obtains 



J«(X,T) = 



1 



1/2 



7rA 2 /5 

/ ds x cxp \{(8 a — 6b)s} exp $ 
■/o 



(-i/fi)^C fc *(iI ) ha C Q 



The explicit expression of <I> ba is shown in (|B4[) , (|B5|) of 
Appendix [El The result is 

2 



$fc« = -- irl V 



& a (T,«)+& o (T,0) 



-^(6 a (T,s) -& o (T,0)) 2 + ^U, 6a (X,T, S ). (50) 



Here Wb a (X, T, s) is given by 



f 2 



w ba (X,T, a) = X(f a - f b )s - f^ s 2 (3T - 2s) 

bm 



/fc2 -{(T-s) 3 -T 3 } 



fafb 



s(T - s) 



6m 2m 
When N — > oo, the integration by s is dominated by the 
stationary point, satisfying 



2m 



(f a -fb)Ts , (2/ a -/ b ) s : 



m 2m 
Taking in advance the limit A — > into consideration 



(51) 



y 6a(r,s) + ^(T,o) _ Q 



(52) 



6a (T, s) - 6a (T, 0) = (/„ - / 6 )(2T - s)s = (53) 

have to be satisfied also. The only solution of (|55|) in the 
range 0<s<Tiss = 0, which also satisfies ([5"T]) . Then 
(|52| becomes X - / b T 2 /2m = 0. Thus, as a function of 
X, $(0)**( 1 ) has a peak at X = &(T). These facts are 
expected; suppose the peak of <3>(°)(X, T) is at In 
order to get non-zero jW(X,T), the peak of #W(X,T) 
shoud also be at £b(T). This can be realized if and only 
if the transition caused by [7Jo] nd from the state \a > 
to \b > occurrs at s = 0. Then the X-integration in 
$(0)*^r(i) i s dominated near X = &(T). To perform the 
calculation, consider the region near s = 0, X — £&COi 

dw ba (X,T, s) 



= 



ds 



(fa -h)[X 



fbT 2 f a — fb 



Ts) + 0(s 2 ). (54) 



2m m 

Thus, near s = 0, the stationry tarjectory s = s(X) and 
the second derivative (which is X-independent) becomes 

m(X-£ b (T)) 



s = s(X) = 



(fa - fb)T 
d 2 UJ ba (X 1 T lS ) (fa-fbfT 



+ 0((X-Zb(T)Y), (55) 



ds 2 



m 



(56) 



Note that the second derivative is a constant. Keeping 
X fixed, s-integrarion is first performed by expanding 
u ba (X, T, s) around s = s(X) 

Wba(X,T,s)=Wba(X,T,s(X)) 

-(T/2m)(f a -f b ) 2 ( S ~ S (X)) 2 + 

In this way, we get 

T 

dsexp(iNuJba(X,T,s)/h) = exp (iXu> ba (X,T,s(X))/h) 



' —2mhm 

N(f a - fb) 2 T 



1 + 0(1/VX) 



(57) 



The above expression is a function of X. Now we con- 
sider its asymptotic functional form when N — > oo. The 
stationary point is given by 







dLU b a(X,T,s(X)) 



(f a -f b )s(X). 



dX 

Therefore X = &(T) and 

oJ ba (X = 6(T),T,s(6(T))) = w&«(6(T),T,0) = 0, 



d 2 CO ba (T, S (X)) 



(fa - fb) 



ds(X) 



dX 2 w " JU ' dX 

Thus for large N, one can write 

exp(iNcj b a(X,T,s(X))/h) 



m 
T' 



exp 



2mhT 
Nm 



6(X-£b(T))+0(l/VN) . (58) 



The factor 1/2 in front appears by the following rea- 
son. By and by s > 0, it follows that X > £ b (T) 
(X < &(T)) if f a > f b (f a < fb) along the station- 
ary trajectory. Therefore, in either case, X = £b(T) is 
the end point of the X-integration and using the formula 
/ °° dx5(x) = 1/2, eq. dSHD follows. 

Other factors in (f5"0|) not multiplied by N are unity for 
large TV, when the stationary value is inserted. In fact, 
consider 

, JiLWffliffl)'} (59) 



x exp 



1 



4A- 



(Zba(T,s(X))-Zba(TM . (60) 



Since ([58]) says that X - &(T) = 0(1/ VN) and s(X) 
X — £b(T), one can estimate for large X as 

Zba(T,s(X))-Zba(T,0) = O(s(X)) 

= 0(X-£b(T)) = 0(l/VN), 

x S(X)) + Cba (T, 0) _ x 

- a a (T, 0)s(X)/2 - &(T, 0)s(X) 2 /4 + • ■ • 
= -C(T,0) S (X) 2 /4+... . 
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This is 0(A-£ b (T)) 2 = 0{1/N). Therefore, both factors 
of ([55]) and (pUj) becomes unity as A — > oo. 

Collecting (|57|) . (f58|) . and adding the term with com- 
plex conjugate, we arrive at 



(x, T) = Y J * b (X, T)¥ b > (X, T) + c.c. 



K. 



(i) 



/7T x 



A 



2ImC 6 *(# ) 6a C 



I /a - /ft 



2m 



0(1/N). 



(61) 
(62) 



This is the first order correction in [i?o] nd to the ideal 
case ((301) • In Appendix [Dl the result (|6"2")l is checked by 
integrating over X first and then by s. Note that y/T in 
([57]) and ([58]) are cancelled, so K b is independent of T 
for each b. In this connection, see Sec lIV Gl 



C. Normalization 

The normalization J dX J(X,T) = 1 leads to 
J dXJ<®(X,T) = 1, 

dXJ {k) (X,T) = (k = 1,2,---)- (63) 



We can check / dXJ^(X,T) = 0. Indedd, note that 

^ ^2hnC*(H ) ba C a 
I 



dxjW(x,T) = y^-Y- 



Here, 1/|/ Q — /b| is real and symmetric under a «-> 6. 
Then, we see that C b (Ho)b a C a /\fa - / b | is a real 

number, so the imaginary part vanishes. 



D. Calculation of j[ 2) (X, T) 

Consider j[ 2 ' of l|2"9")) . It is expressed by 
jf > (X, T) = £ (X, T)*« (X, T) 



1/2 



Y, C *A H 0)a'b{H )baC a ( ds' [ ds 
1 1 Jo Jo 



r/ aa'b 

x expi{-((9 a >s' + 6» 6 (.s - s') - 6» a s)}exp $ a - Q;6 . (64) 
$a'a:b is given in (|C2jl and (fC3|) of Appendix [Cj 



i / u(r,«) + ^ 6 (T, a Q 

<W = " £J ^ 

. ,„ „ m2 iA 



(6a(T, S) - 6a' (T, S')) + — ^aa'; b (T, s), (65) 

(66) 



W a < a;b (T, s) = CJ ba (T, s) - 0J ba ' (T, s') 

= A(Q ba (T, S) - Q ba < (T, s')) + P ba (T, s) - P ba ,(T, s') 



The stationary equation in s is identical to (I51|) : 

(Lw ;b (T,s,s') duj ba (T,s) 
° = * = (6?) 

The solution is written as s = s ba (X). Similarly, we have 

duJ aa '- b (T,S,s') doj ba '(T,s') 







ds' 



'(fa' — fb)[ X — 



ds' 

f b T 2 (f a - - f b )Ts' 



2m 

(2/ a > - h)s' 2 
2m 



(68) 
(69) 



with the solution s' = s' ba ,(X). Now & aa i [b (X, s, s') is 
expanded around s = s ba (X) and s' = s' ba ,(X), and we 
calculate the second derivative at these points, 

M ba = d 2 uj aa '. b /ds 2 = -(f a - f b ) 2 T/m 

~(fa- fb)(Va- fb)s ba (X)/m, (70) 

M ba , = d 2 u aa ,. b /ds' 2 = (f a , - f b ) 2 T/m 

+ (fa' - fb)(2fa' - f b )s' ba ,(X)/m. (71) 

The reult of s, s' integration is 

[ ds [ ds 1 expi{-(6> a ,s' + 8 b (s - s') - 9 a s)} 
Jo Jo 

x exp $ aa / ;6 (X, s,s') 



I (2ir) 2 H 2 
N 2 M ba M ba , 



.exp $ aa , b (X,s ba (X),s ba ,(X)). (72) 



Next task is to study the X-integration. For that pur- 
pose, it is convenient to use the following form for the 
factor appearing in (|65[k 



exp 



1 f Y Zba(T,s)+W(T,s') 
M 2 



x exp 



1 



4A : 



■(^ b a(T,s)-^ ba '(T,s')f 



exp 



2A 2 \ X ~^a(T, S) 



x exp 



2A 2 



X-Z ba ,(T,s') 



(73) 



Due to the structure of & aa <;b(X, s ba (X), s' ba ,(X)), the 
resulting dependence on X differs for a' = a and a' ^ a. 

The case a = a' 

uj aa . b (X,s ba (X),s' ba ,(X)) = holds since for a = a' 
s ba (X) = s' ba ,(X) = s(X). Thus the result of the X- 
integration is a constant independent of A. Consider 
the factor (|T3"]) contained in $ a ^.j of (f6"5)) . Inserting the 
stationary value s(X) = s ba (X) in the first factor of the 
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right-hand side of (|73)l . we concentrate on X — £b a (s(X)). 
The factor in the second parenthesis in (|5ip is rearranged 
as 



X 



1 

2m 



f b T 2 + 2(f a - f b )T S (X) - (2f a - f b )s(X'f 



1 



= X-i ba {T,s{X)) + --f aS {Xf. 

2m 



In this way, we get 

X-i ba {T,s{X)) 



-{f a /2m)s{Xf 



(74) 



(75) 



Since we are considering A — > (after N — > oo), X — 
£ 6o (T,s(X)) -> 0, implying s(X) -> 0. On the other 
hand, by ([4"6")l . one can approximate 



X-& a (T )S (X)) 



= X - 



fbT 2 (f a - f b )T 



s{X) = 0. 



2m m 

Soving this relation for s(X) and inserting it back into 
([751. we conclude 



X-& a (T, S (X)) = 
Thus we obtain 

P = cxp 



/a 



2m (/„ - / b ) 2 T 2 



X 



fbT 2 
2m 



bo 



Cba = 



A 2 T 4 
ifl 



X - 



2m 



4m 2 (/ a -/b) 4 ' 
Here the following formula is adopted. With C > 0, 



lim exp 



-C 



^S(X-a), 



where 7 = J™ exp(— z 4 ). In this way, we arrive at 

fbT 2 ' 



2m 



The factor 1 /2 is present for the same reason as given con- 
cerning ([58)1 . Collecting all factors, the result for J^ a=a , 
is obtained as follows. In doing so, Sb a (X) and Sb a '(X) 
appearing in M bQ , Mb a > of ([70jl . ([71)1 can be set to zero, 
which is inserted into ([72jl . We use ([37)) and set p = 1. 

</S=a< = E ^(-i/ft)(i/S)C:,( J ff )a' 6 (^0)6aC 
a— a'6 

TVm 



x 



1 



2tt 



ttA 2 2-kKF y/JjTF 
1 Trh2m 



N{f a -fb) 2 T cT 2 



^ T1 s(x- fbT2 



2m 



4i> 



(2); 



/ b r 2 

2m 



K. 



(2) 



V2 7% M^ y, C*(F )a6(-ffo)6aC 



l/a - /&IV|/a| 



Aft 



(76) 



The case / a' 



In this case, s'(X) 7^ s(Jf) so $ a >a;b(X, s(X),s'(X)) does 
not vanish and is a function of X of O(N), producing a 
stationary phase. We need to pin down the position. 
(We write s ba (X) = s(X) and s' bab (X) = s'(X).) By the 
stationarity in s and s', eq.(J75J and a similar equation 
for s' with a replaced by a' hold; 



X-Sba'(T,s'(X)) 



f 1T 2 

2m v ; 



(77) 



Eqs. ([75"j) and ([77]) assure that we can limit our discussions 
near s(X) = and s'(X) = 0, as stated just below ([75)1 . 

Now, using the stationarity in s and s', the stationary 
condition of $ a 'a;b(X, s(X), s'(X)) with respect to X can 
be written as, 

d* a , a]b (X,s(X),s'(X)) _ d* a , a]b (X,s(X), S '(X)) 
dX dX 
= (f a s(X) + f b {T - s{X)) - ((f a ,s'(X) + f b (T - s'(X))) 
= (f a -fb) s (X)-(f al -fb)s'(X). 

Higher derivatives are obtained by differntiating the sta- 
tionary equation of s or s' . Differentiate ([5T)1 by X; 

1 d 2 s 
-2— (/„-/.,)— T 

Setting s = in the first eqution, (f a — f b )(ds/dX) — 
m/T, which is inserted into the second. Thus one gets 



{fa - fb) 



d 2 s m 2 (2/„ - f b ) 



dX* T3 (/„ - f b ) 2 



Using the similar equations for s', with the replacement 
a a', we get at s(X) = 0, s'(X) = 0, 



d$a>a;b 



dX 

1 rf 3< £q'a; 

3! dJf 3 



0. 



dX 2 



0. 



1 d 2 s d 2 s' 

3! (/ ° " /fc) rfX 2 " ~ (/a ' ~ ^dX 2 



m 2 {fa — fa')(fa.fb + fa'fb — Zfafa') _ C a > a b 



6T3 (fa-fb) 2 (fa>-fb) 2 

In this way, $ a 'a:b becomes 



T 3 



<$>a'a;b 



iN C a , ah ( x f b T 



2\ 3 



h T 3 \" 2m 
Here we use the formula of symmetric integration 
lim exp iaN(X — a) 3 = lim cos aN(X — a) z 



(78) 



(79) 



\a\N)- 1/3 j'S(X -a), 7' 



dz cos z. 



10 



Just as in the discussions for given in (|59|) and (|60)l . 
the factors which are independent of N become unity, 
since N — > oo is taken before A — > 0. The result for 
a =/= a' is thus obtained; 

j (2) - V 

a^a' \b' 



{-i/h)(i/h)C*,(Ho)a>b(Ho)baC a 



1 Nm 2ir 



ttA 2 2irHT y/DD* 
1 7rft2m 



N\f*~ h\\f a > ~ h\T CV*WI* 



2m 



NN 1 / 



tfi/aZ^b d l A 2m 



A", 



/(2) _ VmVi ^ p ^ C;,(ff ) a , fc (ff ) to C' o fi 1 ^ (8Q) 



Aft 



\fa - fb\\fa> ~ fb\C a , ab 



Here Re signifies the real part. Under the exchange a «-> 

1/3 

a', both C*,(H )a'b{H )baC a and C a / ab become complex 
conjugate. So the summation over all a ^ a 1 is equivalent 
to apply Re. The results (|76|) and (|80|) have been checked 
by interating by X first and then by s, s'. 



E. Calculation of 



The remaining term of 0(([Ho] nd ) 2 ) appearing in 
is ^°>(b,X,T)^ 2 \b,X,T)+c.c, see ®. Actually, it 
is not necessary to calculate this term from the start. 
One can invoke to the normalization condition (16311 with 



(2) (2) 

k = 2. Consider the result of J{ . If we write J{ = 
E 6 jf )& , then 

J {2)b (X,T) = * (1) *(&,A,T)<I< (1) (&,A,T) 

= J2 Vaia;bC* a ,(H )a>b{Ho)baC a 8(X - &(T)). (81) 

a, a' 

V^' a; 6 is given by the sum of ([To) and ([50)) and has dif- 
ferent forms for a — a' and a ^= a . As dictated from the 
definition of of (23J), we can write J^ 2) = £ 6 4 2)b , 
where J^ 2 ^ has the form 

jf )& (A,T) = ¥°>(b,X,T)^ (2 \b,X,T) + c.c. 

= Wba-cC* b (Ho)bc(Ho)caC a S(X - &(T)). (82) 

c,a 

Note the differnce of the index structure of b between 
flUJ and ([82"]). Inserting these two into {63]) (fc = 2), 

Wba-cC b (Ho)bc(Ho)caC a 

c,a,b 

= — ^ V a >a;bC*,(Ho)a'b(Ho)baC a 
a 7 a' ,b 

is obtained. After renaming the index of W, and recall- 
ing that the above equation holds for any C a and any 



operator Ho, W a ' a ;b = — V a ', a ;b follows. Thus we arrive 
at for a = a' and a ^ a 1 separately; 



t(2) (y T\- 1 ^7V^F 

J 2;a=a'\ X i 1 )- — 



E 

a , 6 



AT VAft 

C* (i/o )ab (Ho)ba C a 



\fa ~ fb\y/Ja 

1 2^'mypK 



S(X-UT)), 



(83) 



AWV3 Aft 



(84) 



a^a' ,b 



I /a' /fell /a. /ft I C a ' a6 



F. The result for J(X,T) up to ([#o] nd ) 2 

Collecting all the results of gDJ), fH, f7B jl .l|5D )l fSJ) 
and (JSU), the signal function up to 0([H o ] ad ) 2 is 

J(X,T) = J(°)(A,T) + J^(X,T) + j[ 2) a=a ,(X,T) 

+4%a> (X, T) + J% =a , (A, T) + 4% a , (A, T) 



E 

6 



IGJ 2 



Ab 
N 



B h 



NN 1 / 3 



6{X - f b T 2 /2m). (85) 



Here each coefficient is given by 



A bl + A b2 , 

^ 21m C^(H ) ba C a 
A ^ 



1,2 



V2jy 



Aft 



I /a - fb\ 

{Ho ) afc (Hq )baC a 



(86) 



E 



C b (Ho)cb(Ho)ca,Cb 



I /h /c I v/o 



(87) 



Si, 



2^)' rriypKh 



1/3 



Aft 



x Re 



Ca,(Ho)a'b(Ho)baC a 



a^a' I /a — /b| l/a 



C b (Ho)bc (Ho)caCa 



b^a,c 



\fb 



/c||/a fc\C bac 



(88) 



The origin of the power of 1/N is the stationary phase 
integrations by s, s' and A. Each s-integration brings us 
1/\N, but the results of the A-integration (or equiva- 
lently of the formula ©) depend on the situation; 1/yN 
for Abi, (1/VN)° for A b 2 (no stationary phase) and for 
B b , 1/A 1 / 3 . In Appendix [El the order of neglected terms 
are estimated. They are shown to be down by at least 
1/y/N compared with (|85|) . 

As T -> 0, above J(A, T) goes over to J(A, 0) = S(X) 
by the same relation that follows from the normalization 
condition. Note further that when [Ho,0] = 0, all the 
correction terms are absent and the ideal result becomes 
exact up to the order considered here. 
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G. T-independence of Abi,b2 and Bb 

Above correction terms A b i, A b 2 and B b are inde- 
pendent of T, besides \C b \ 2 - The reason is simple; as 
we have seen above, the fluctuation into any channel 
b occurs at t = for N — •> oo, and after that the 
state develops by ?7jv,6 01 (H3) and ifTS]). so the claasical 
time evolution is realized. Consider the signal function 
J(b,X,T) = \$(b,X,T)\ 2 defined for each channel, i.e. 
one term in the sum JSHJ of J{X,T) = J2 b J(b,X,T). 
Then once the above fluctuation at t — is taken into 
account, J(b,X,T) evolves as 8{X — £b(T)) without any 
further T dependence. Therefore, when we integrate 
J(b, X, T) in R b of (gTJ and define 

J(b,T)= I J(b,X,T) = I J(X,T)dX, 

then J(b, T) is independent of T. Since this holds for any 
Ho or for any parameters in the thory, we conclude that 
Mli M2 or B b does not depend on T. This statement 
applies even if one includes higher order of Hq- However, 
when the diffusion processe is taken into account, the 
above assertion does not hold. 



larger than ([89)) and may be the main correction to the 
ideal zeroth order result of ([55)) . The reason why this 
term is large is, as stated before, understood if we look 
at dnSJ). When a = a', u> aa;b (X,s ab (X),s' ab (X)) = 0; 
after inserting the stationary condition of s and s', the 
exponent of <& becomes O(l), not O(N), so the cancela- 
tion of the phase as a function of X does not occur and 
the resulting value becomes large. When, for example, 
N = 10 5 particles are participating in giving the signal, 
the corrction of the order 10~ 7 is expected. This can be 
within the experimental confirmation. 

Finally as for B b , we note Cbac/h ~ m2 1 fa by ([75)1 . In 
this way, we get 

B b 1 ma 2 1 a 

awV3 = jvjvvs ~Y~ x (c o6c /fi) 1/3 x A 

= 1 ( faa \ - ~ 1 x 2 5 x 1(T 4 
NNVsytP/ima 2 )] A NN 1 / 3 

This is a small number. 



V. DISCUSSIONS 



H. Numerical estimates 

Here we try to estaimate the order of numerical values 
of the correction terms. Consider first A b i given in ([55)k 
The numerator is the order of enegy of the object. As dis- 
cussed in deriving (|38|) . f a a order of the energy a partiacle 
in the detector recerives from the object. Thus we regard 
C b *(H )baC a is is of the same oredr as |/ - f b \a - f a a. 
In this way, adopting A = 10 _3 cm, a = 10~ 7 cm, 



Mi }_ a_ 
~N~ ~ N X A 



1 

N 



x 10" 



(89) 



The condition a/NA « 1 is our requirement for the 
measuring process to be sensible. Thus one has to require 
N ;» 10 4 . We have stressed thoughout the paper that 
the limit N — > oo is kaken before A — > 0, which is denoted 
by =>. The requirement a/NA <<C 1 agrees with this limit 
and at the same time gives the precise condition of this 
limit. In our case, the "macroscopic" number of particles 
in the detector required to produce experimental signal is 
expeessed as N ^>> 10 4 , which can be much smaller than, 
for example, the Avogadro number. The last conclusion 
is stated frequently for realistic devices. 

Next, we estimate A b 2 of ([57)) . Similarly as above, the 
factor C b *{Ho)cb{Ho) C aCb is regarded as the same order 
with |/ Q — fblfaO 2 ~ fa - 2 - Then apart from the numerical 
constant, A b2 ~ (a 2 /K)^\f a \m/ A. In this way, 



— x 4.0 x 10^.(90) 




In order to see whether or not our result ([85)1 is specific 
to our model, let us consider the general Hamiltonian (JSJ, 
combined with ^ and apply the arguments of Appendix 
|A1 When N —> oo, the leading term agrees with the ideal 
result of ([55)) . This is because Ho is negligible and the 
object stays in the prescribed eigen-state \b > of O. For 
fixed 6, the detector evolves by Hd + Hj(Xb, [xi]), so the 
stationary path X^*(t) depends on b. By normalization, 
the signal function ()A2p applied at t takes the form 



Here we have used the same values as in (1381) . For the 
electron, h 2 /(ma 2 ) ~ 10~ 13 erg. The above l[9"0")) is much 



J(X,i)=^|C 6 | 2 5(X-Xf( i )) 

b 

for general case also. Expanding in [Ho] nd , we use (|14|) . 
with one s integration accompanying each Ho- Since 
each s integral is dominated by the stationary phase, 
one s brings us 1/%/A, irrespective of the form of h 
and g. We conclude that k-th order correction terms 
containts at least the factor [\/^/N) k . Whether extra 
factors of X/N appears or not depends on the form of 
H D + Hj(X b , [xi]), just as in our model extra factors of 
1/N appeared through A-interation. This was pointed 
out just below ([88)) . We have to study more realistic 
detection process and fix the above correction terms, in- 
cluding the numerical estimates. In doing so, some sim- 
plification of the Hamiltonian will be required of course 
to make the problem tractable. 

Our scheme can also be applied directly to the system 
described by the quantized field. Indeed, the field theory 
is much more suited to handle the macroscopic system, 
especially when taking the thermodynamic limit. Also, 
the measurement theory in the relativistic case can be 
studied using the field theory, since the formalism of the 
relativistic field theory is firmly established. 



12 



Finally, the most difficult problem of reduction, or ein- 
selection, is left out of the discussions in this paper. The 
extension of the dynamical reduction thory [l5[ to the 
relativistic case is controversial. However, in our case 
the application of the stationary phase to the relativistic 
field theory is straightforward. Also our results in this 
paper are independent of the precise mechanism of the 
reduction process, since we rely solely on the Schrodinger 
equation, which holds independently of how the reduc- 
tion is realized. 



APPENDIX A: NON-SEPARABLE CASE 

To obtain the equation of motion for a macrovariablc 
is the same problem of how to get the effective thory of 
a collective mode in many particle system. There are 
several methods but here we select the one which con- 
sists of inserting a delta-function in the path-integral. In 
the limit N — > oo, it becomes equal to the method of 
the Legendre transformation [22T |. We consider the case 
where X is given in collective notation [xi] as g([xi])/N. 
Extension to the case where g includes [pi] is not difficult. 

The signal function 

Let us write the wave function of the detector plus object 
system at some fixed time as x). For any thermo- 

dynamically normal macroscopic system, the wave func- 
tion has the form [2111 



*([Xi 



Ge 



(Al) 



where F = F(\xi] , x) is an extensive quantity and is of the 
order O(N), so we write it as F([xi],x) — NT{\xi\, x). 
On the other hand, G — G([xi],x) describes the mi- 
croscopic details. By normalizability, the real part of 
J-([xi], x) is positive definite and when ^"([xj], a;) changes 
by a finite amount, it describes macroscopically different 
state of the system. When we sum up wave functions hav- 
ing different values of T, it represents the mixed state. 
The pure state is obtained by summing up various G"s, 
with the same T . The signal function defined in ([25]) can 
be generalized to the non-separable case as follows. 

J{X) = N J dx /" JJ dxi 

i 

xmixilx^SiNX-gdxi])) (A2) 
= N J dx j Y\ dx i J dj 

x exp f - N2BeJ r {[x i ], x) + i{NX - g{[xi])}j 

/OO 
dj exp(-JTO') + iNjX) = e~ K( - x \ (A3) 
-oo 

We have dropped microscopic G term since it does not 
affect the stationary condition for X. Note here that 
J(X) is real by definition, so K(X) is a real quantity 



and positive definite. Now H(j) is extensive, so the func- 
tion H(j) and K (X) are both extensive proportional to 
N. The results (JXTJ and |M| have been obtained HJ 
by adoping the field theory to describe the macroscopic 
system, which automatically takes into account the sta- 
tistical factor l/N\. The essential condition for these two 
equations to hold is the short range character of the inter- 
ation among constituent particles in the macroscopic sys- 
tem. Although the field theoretical approach is not taken 
in this paper, here we apply (|Alj) and (|A3j) to the quan- 
tum mechanical iV-particle system, since these results 
express the general property of the thermodynamically 
normal system. Writing K(X) = NK(X), we first ex- 
pand fC(X) around the the stationarity solution X = X st 
satisfying dJC(X) / dX = 0. Suppose the second deriva- 
tive at the stationary point is positive; K^ 2 \X st ) > 0. 
(Since IC(X) is positive definite, at least one stationary 
point with > exists. The solution with /C (2 ) < 
represents an unstable state.) In the limit N — > oo, the 
analog of © for the case of positive definite tC(X) holds; 



lim e-^W 



-NK(X" 



2tt 



N- 



-NK."(X S 



-5(X ~X st ). 



This is equal to S(X — X st ) by normalization. 

When J(X) has the distribution in A, it represents 
the mixed state. Indeed, the signal function is rewritten 
as 



J{X) = / dY J(Y)S(X - Y) 



dYJ(Y) v / N/?:exp(-N{X - Y) 2 ). 



lim 



In this form, J(X) of ([12]) is seen to be a superposi- 
tion of different JF's, with the weight J(Y). In real 
detection process, the detector is so arranged that T 
depends on the eigen-value A a of the object opera- 
tor O to be measured. Writing the time t explicitly, 
F{X a ,[x t ],t) = NF(\a,[xi],t) and K = NIC in (TM| 
becomes NK.(X a ,X,t). Then the stationary solution of 
dIC/dX = depends on t and a; X = Xjf (*). For the 
model of (JTSJ), this is written as £ a (T) in (152]) . 

The equation of motion of A st (t) itself can be obtained 
by the method of double paths Legendre transform, as 
discussed in detail in 1 22] . 



APPENDIX B: CALCULTION OF $ 6a 

By (EH), (H]) and (08]), jW(X,T) is given by 

(A, T) = Y J y {0) *(b, X, T)#W (6, A, T) + c.c. 



7tA 2 /9 
T 



(-i/h)J2c*b(Ho) ba C a 

ab 

/ ds x exp i{(0 a - 9 b )s} exp <I>f, a + c.c. (Bl) 
Jo 
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By X - &(T) = R ba (X, T, 0), $ 6q is expressed by 
N 2 m 2 fR ba (X,T,s) 2 R ba (X,T,0) 2 



$ba. = - 



2h 2 T 2 \ D 



D* 



Nmi 
i 



{R ba (X,T,s) 2 -R ba (X,T,0) 2 } 



N{XQ ba (T,s) + P ba (T,s) (B2) 
11 

-XQ ba (T,0)-P ba (T,0)}. 
Here we use (I3T1) and rewrite the first term of $> ba as 
N 2 m 2 fR ba (X,T,s) 2 R ba (X,T,0) 2 ' 



2h 2 T 2 



D 



D* 



pA 2 
1 

iNm 



(e 6a (r, S )-e fe a(T,o)) 2 



(B3) 



" 6a(T ' s))2 ~ (x ~ 6a(T ' 0))2} - 

Approximating p — 1 , (f>f, a is now rewritten as 



^ Q (T, S )+6a(T, 0) 



- 4^2 (&» ( T > s ) ~ 6a (T, 0) ) 2 + ^ ^ ba (X, T, s), (B4) 
with to ba (X, T, s) defined by 

u ba (X, T, s) = X(Q(T, s) - Q(T, 0)) + P(T, s) - P(T, 0) 



X(f a - f b )s--^-s 2 (3T-2s) 

fb {(T _ s) 3 _ T 3 } _ fah s{T _ s) 2 



6m 



2m 



(B5) 



APPENDIX C: CALCULATION OF $„,„ 



By (|48|) . the explicit form of <& a 'a;b defined in (jM)) is 
N 2 m 2 fR ba (T,s) 2 , i? b ^(T, S ') 2S 



$a'a;b 



2ft 2 T 2 V D 



Nmi 
2HT 



{Rba(T, 3 y -R ba >{T, s'Y) 



+ -N{XQ ab (T,s)+P ab (T,s) 
n 

-XQ a , b (T,s')-P a , b (T,s')}. (CI) 
As in (|B3|) . we rewrite 

N 2 m 2 ^ 2 ofr s\2* 



2h 2 T 2 



( R(T,s) 2 R(T,s') 2 \ 
\ D D* ) 



pA 2 



+ i ^(6a(T, S )-6a'(r,.s')) 2 



iNm 
2phT 



{(i-ur,^-(i-^(T, s ')n. 



By setting p = 1, $> a 'a;b is written as 

ffi 1 L 6a(J»+6a'(JV) 

-^2 (6a(T, s) - Cba'(T, s')) 2 + ^w a /«;6(T, a), (C2) 

Ua>a;b{T, S) = LO ba (T , s) - LO ba < (T, s') (C3) 

= X(Q 6a (T, s) - Q bo ,(T, s')) + s) - P ba ,(T, s'). 



APPENDIX D: INTEGRATION BY X 



Since each term in J2 b °f (ED is given by K^S(X — 
£b(T)), is calculable by integrating over X first and 
then by s for fixed &. After integration by X, the term 



exp 



-1 j Y S ba (T,s)+Z ba (T,0) Y iN 

A>r 2 — } +T x(fa - fb)s 



changes into 



y / nA 2 p 



exp 



A 2 A 2 p 
Ah 2 



2„2 



(/a - fbYs 



Ni 



+ ^ ( 6a (T, S) + Cba(T, 0) ) (/„ - f b )s 



Then the dominant region of s-integration is s = 0(X/N). 

The second term in [ ] can be shown to be 0(s 3 ), if 

it is combined with P ba (T,s) — P ba (T, 0). Indeed, 

(&a(T, s) + 6a (T, 0))(/ o -f b )s + P ba (T, s) - P ba (T, 0) 
~ (fa- fb)(2fa- fb)s 3 /6m. 

for small s. By TVs = y, s-integration now becomes 

1 fTN 

— J dyxcxpi{(6 a -9 b )y/N}^^ 
A 2 P/ 



x exp 



(fa - fb)V 



Ah 2 



0(l/N 2 ). 



It gives hir/(N\f a - M)(l + OQ./N)). Thus we get 

-iV^E a C* b (Ho) ba C a 



NA\f a -f b \ 
Adding term with c.c, we get (|6"Tj) and (|G*2"]) . 

APPENDIX E: NEGLECTED TERMS FOR JV -> oo 

First, let us restrict the arguments up to 0([Ho] nd ) 2 
and consider the correction to |C b | 2 , A b and B b of ([55]) . 
1. |Ch| 2 ; As stated in Sec lIII El the correction comes 
from the diffusion. Including the diffusion, A is re- 
placed by yfpA. Since <fp = 1 + 0(1/A 2 ), |C 6 | 2 changes 
into | Cfo | 2 ( 1 + 0(l/N 2 )). Thus the neglected terms are 
0(1/N 2 ). The fact that it is small numerically has been 
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checked in Sec lIIIEl 

2. A b \ The correction comes from the higher order of the 
fluctuation around the stationary path, which is down 
by 0(1/ >/~N) compared with the term retained. In ob- 
taining Abi of (f86|) . both s and X have fluctuations of 
0(1/ VN). As for A2b, s and s' has the same size of the 
fluctuations but X has no stationary point. Therefore 
Abi.2 are replaced by A b i^(l + 0(1/Vn). The leading 
order of neglected terms in Ab/N is thus 0(1/N^/N). 

3. Bb\ Besides s and s' having the fluctuation of 
0(1 /yN), X fluctuates near the stationary phase with 
the size of Oil/N 1 / 3 ), see (JTHJ). Thus B b is replaced by 



B b (l + Oil/N 1 / 3 ) and the corrections to Bb/NN 1 / 3 is 

(9(l/7ViV 1 /3)(l/iV 1 /3). 

Next, consider the term ([Ho] nd ) k for k > 3. In the 
expansion (H3J), one factor of [i/o] nd accompanies an in- 
tegration over the parameter s. Integration around the 
stationary phase produces one factor 1/yN, so (l/y/N) k 
appears. Extra factor of some power of 1/./V may appear 
depending on the result of the X-intergation, as stated 
above. Therefore, we can say that ([Ho) nd ) k term has 
the power of at least (l/VN) k . 

Summarizeibg, eq. (|85p is correct up to the order re- 
tained there, the leading correction being 0(\/N^/N). 
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